Charmed-meson spectroscopy in QCD sum rule 
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We elaborate masses of open-charm mesons (cd and cs) with J F = 0~,1~,0 + ,1 + from a view- 
point of ordinary light-heavy systems in the analysis of standard Borel-transformed QCD sum rule. 
This analysis is implemented with the operator product expansion up to dimension 6, permitting 
corrections to the order q s and to the order m s , and without relying on l/m c -expansion. The 
obtained results following our stringent criteria for the continuum-threshold dependence, the Borel 
window and the Borel stability, indicate that the charmed-meson masses in the + channel are 
overestimated by 100 ~ 200 MeV in comparison with the experimental data, which were recently 
reported as the rather broad state of D* + (2351) and the extremely narrow state of D* + (2317), 
respectively. Such large mass-discrepancies from the data cannot be seen in other channels, where 
conversely our results of the cs-meson masses are even underestimated somewhat in comparison 
with data, independent of the value of the strange-quark mass adopted in our calculations. From 
these results, it might be expected that the measured low mass of D| + (2317) is a manifestation of 
an exotic state with the structures of a four-quark or a molecule, while at present the D* + is not in 
conflict with the existing data due to its large width. 

PACS numbers: 14.40.Lb, 12.38.-t, 11.55.Hx, 11.55.Fv 



I. INTRODUCTION 



The detailed study of charmed-meson spectroscopy has 
provided us with understanding of not only its struc- 
ture and property peculiar to the light-heavy system, 
but also widely the basic property of the strong inter- 
actions. Especially, the charmed mesons have played a 
part for the success of a potential model analysis in meson 
spectroscopy through a unified framework However, 
recent experimental discovery of scalar and axial-vector 
cs-mesons 0, SJjl HI El i which have been only missing 
p-wave states Q, cast a shadow on such a success of 
the potential models, because their predictions [cl lil IT(i| 
give heavier masses by 100 ~ 200 MeV than the obser- 
vations. Similar results have been seen also in QCD sum 
rule coupled with heavy quark effective theory (HQET) 
[Tl| and lattice QCD simulations fl2l . In the mean- 
while, an interpretation of states with J p = + ,l + as 
the chiral partner of ground states with 0~,1 ~~ in the 
heavy quark limit requires relatively lower masses for the 
(0 + , 1 + ) states than the above predictions. The analysis 
based on such idea gives the results comparable to the 
observations [l4l llj| . This problem with the mass dis- 
crepancy has triggered a lot of debates on the natures 
of the cs mesons in the viewpoints of modification of 
theoretical models, e.g. involving Z^iv-mixing lfj, or in- 
troduction of unconventional structures like two-meson 
molecular states ^3 113 or four-quark states [Til l2fj| or 
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their mixture |2lJ. There is, however, still no definite 
explanation for the discrepancy. 

In this paper, we study eight states of the charmed 
mesons with J p = 0~,1~,0 + ,1 + , which are composed 
of the corresponding four states of cs-mesons and simi- 
larly four states of cn-mesons with n = u, d. The masses 
of their states are elaborated through all the same analy- 
sis from the viewpoint of the conventional interpretation 
as light-heavy-quark bound systems. We then focus at- 
tention on the masses of two cs-mesons with + and 1 + 
in comparison with the experimental data. This method 
would be quite useful to extract information on specific 
states like + and 1 + in the analysis involving relatively 
large theoretical uncertainties. Here, in the 1 + channel, 
we deal with only a channel of charge- conjugation par- 
ity + 1, i.e. J p ( c ^ = channel, although indeed the 
light-heavy system does not have this parity as a good 
quantum number and thus takes place a mixing between 
l+v+) and channels. 

In the present work, such a physical motivation can be 
realized technically by relying on the standard QCD sum 
rule (QSR) approach |22l |23 |. which evaluates hadron 
properties by using the correlator of the quark currents 
over the physical vacuum. We perform its calculation 
without taking l/m c -expansion, where m c is the charm 
quark mass. It is because the charm quark is compara- 
tively light and some properties of charmed-meson sys- 
tems seem to be not necessarily valid for such an expan- 
sion, especially for their masses. However, some analyses 
have approximated at the first few terms of the expan- 
sion based on the HQET U [H [12 . 0ur anal y sis is 
implemented with the operator product expansion up to 
dimension 6, permitting corrections to the order a s and 
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to the order m s , where a s is strong coupling constant and 
m s the strange-quark mass. In Borel-transformcd QSR, 
so-called Borel sum rule (BSR) j^, we extract charmed- 
meson masses with imposing criteria for the continuum- 
threshold dependence, the Borel window and the Borel 
stability. Indeed, we find that the D(0+), D s (0+) -meson 
masses are overestimated by 100 ~ 200 MeV compared 
with the experimental data, while such a tendency cannot 
be seen in other channels. Conversely, in other channels 
the cs-meson masses are even underestimated somewhat 
in comparison with the data, independent of the value 
of the strange-quark mass adopted in our analysis. Ac- 
counting for the measured broad width of the D(0 + )- 
meson, the state seems to be not in conflict with the ex- 
isting data. From these results, it is quite possible that 
the measured low mass of D* + (2317) is a manifestation 
of any exotic states with the structures of a four-quark 
or a molecule. 

This paper is organized as follows. In sec. II, we 
first mention the current status of charmed-meson spec- 
troscopy in more detail, from both experimental and the- 
oretical aspects. In sec. Ill, we formulate the general 
form of BSR and derive the analytic forms of cs-mesons 
in four channels of 0~, 1~,0 + , 1 + , respectively, and sim- 
ilar relations for the cn mesons are also easily obtained. 
Our criteria requested in the analysis of BSR are given 
in sec. IV and numerical results obtained thus are also 
shown. Finally, summary and discussion are devoted in 
sec. V. 



II. CHARMED-MESON SPECTROSCOPY 

In general, the classification of mesons containing a 
single heavy quark is interpreted with the help of heavy- 
quark symmetry (HQS) 7], i.e. the symmetry valid 
for the infinitely heavy mass of charm quark. Under 
this symmetry, the strong interaction conserves total an- 
gular momentum of the light quark, j. In the mean- 
while, total angular momentum of the light-heavy sys- 
tem, J, should be still regarded as a good quantum num- 
ber of the system, even if the HQS breaks down. In- 
deed, the charm quark is much heavier than the QCD 
scale (Aqcd — 0.25 GeV) but not infinitely heavy. In 
this way, we explain the classification of charmed mesons 
in terms of (L, S, J, j) appropriately according to cir- 
cumstances. Here L and S denote the orbital angu- 
lar momentum between the light and heavy quarks and 
total spin of the system, respectively. For instance, 
two ground states (L = 0, J p(c) = 0- (+) ,l~ ( ^) form 
the j = 1/2 doublets and four the first excited states 
(L = 1, J p ( c ) = 0+(+\l+(+\l+(-),2+(+)) can be un- 
derstood as j = 1/2 doublets (J p ^ = 0+(+\ l+M) and 
j = 3/2 doublets (J p ^ = ]+(-), 2+(+>). The states 
up to L = 1 with two classifications are shown in the 
horizontal axis of Fig. 2] 

Four of these six states, i.e. j = 1/2 doublets at L = 
and j — 3/2 doublets at L = 1, have been observed 



over the past two decades f~ 1975 — 1994) following the 
discovery of open charms H, because their states have 
relatively narrow widths |2J|. The properties of masses 
and widths for both charmed non-strange and strange 
mesons match well with thepredictions from the poten- 
tial model analyses 0, H Ej. Recently (2000 ~), the 
leading candidates for two missing states, j = 1/2 dou- 
blets at L = 1, have been also discovered: the first ob- 
servation of D® with the features of J p = 1 + was re- 
ported by the CLEO Collaboration |25| and it was also 
observed by the BELLE Collaboration^^ ■ The D\ has 
the mass M D o ~ 2423 MeV and the very broad width 
r^o ~ 329 MeV, where we adopted the average of the 
above two data. Besides the D®, the BELLE Collabora- 
tion observed also a very broad resonance Dq° with the 
features of J p = + [2g]. Though a marginally com- 
patible result with that of the BELLE for the mass, the 
FOCUS Collaboration also confirmed the evidence of the 
Dl° and at the same time observed its charged partner 
Dq + , which also has the rather broad width [27] . The av- 
erage of the BELLE and FOCUS data on the D* Q ° gives 
M n ,o ~ 2351 MeV and r D .o ~ 262 MeV. The FOCUS 
data on the DX + is M n ,+ ~ 2403 MeV and r n ,+ ~ 283 
MeV. As seen from these data, the j = 1/2 doublets at 
L = 1 have the very broad widths of several hundreds 
MeV in contrast with the j = 3/2 doublets at L = 1. 
Their masses and broad widths are not in conflict with 
predictions from the potential model analyses 8, 9] . Very 
recently, also for the charmed-strange mesons four ex- 
perimental groups (BABAR, CLEO, BELLE, FOCUS) 
independently have observed the leading candidates for 
missing j = 1/2 doublets at L — 1. Contrary to our natu- 
ral speculation about their masses and widths in analogy 
with the charmed non-strange mesons, surprising facts 
were reported by the observations of their groups, as 
mentioned below in detail. 

The £)*^"(2317) was first discovered in the Dfn° de- 
cay channel with the very narrow width by the BABAR 
Collaboration 0], and in the same decay channel its ex- 
istence was soon later confirmed by the CLEO Q and 
the BELLE 5] Collaborations. Its mass was measured 
to be 2317 MeV commonly, which is approximately 46 
MeV below the DK threshold. The width was deter- 
mined to be less than 4.6 MeV, which is compatible with 
the experimental resolution. Here we adopted the sever- 
est upper bound from the BELLE |5| . More recently, the 
FOCUS Collaboration reported a preliminary result 
on the measurement of the mass 2323 MeV, slightly larger 
than those from the other three experiments. Together 
with the Z)*+(2317), the CLEO and the BELLE [| 
observed another narrow resonance in the D*+tt° decay 
channel, with the mass close to 2.46 GeV, where we de- 
note the state by D'X- Later, the BABAR confirmed its 
existence in the same decay channel 3], although they 
have already found such a signal in Ref. [2| ■ It can be as- 
signed to J p = 1 + from the observation of the radiative 
decay in Dfj and the analysis of helicity distributions. 
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Their average mass is 2459 MeV, which is approximately 
46 MeV below D*+K° and the width is less than 5.5 
MeV, which is the severest upper bound observed by the 
BELLE 0. 

On the other hand, prior to such observations, most 
potential models as typified by Refs. have predicted 
that two missing states of + and 1+ (L = l,j = 1/2 dou- 
blets) should be massive enough to decay into DK and 
D* K in a s-wave, respectively. The widths are then ex- 
pected to be very broad due to the strong decays. The 
spectroscopy of cn and cs mesons with L = 0, 1 is drawn 
in Fig.^in order to compare theoretical predictions (solid 
lines for the cs and dashed lines for the cn) with the ex- 
perimental data (closed circle for the cs and open circle 
for the cn). We refer to numerical results in Ref. || as 
a typical calculation by the potential models. Their val- 
ues on the masses are written down in the figure. After 
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FIG. 1: Spectroscopy of cn and cs mesons with L — 0, 1. 
The solid lines for the cs and dashed lines for the cn show 
numerical results quoted from Ref. @ and the experimental 
data to be compared are given by closed and open circles 
HSQBHIslliilglllia, respectively. Their values on 
the masses are written down in this figure. For the purpose 
of reference, two dash-dotted lines display the thresholds of 
D+K° and D*+K° decays for the cs mesons. 

the observations, using data of their new states, Cahn 
and Jackson 0] showed that there are no model param- 
eters to explain all the masses and widths of the + and 
1 + states consistently within the potential models. Thus, 
large discrepancies between the observations and the the- 
oretical predictions for the masses and the widths have 



renewed theoretical interests in the issue of the charmed- 
strange meson spectroscopy, especially for the + state 
|28j . Such a large discrepancy of the mass in the 0+ 
channel is seen also in the QCD sum rule analysis based 
on HQET |Ti) and lattice QCD calculations in the static 
limit ^3 or with dynamical quarks [l3|. although their 
analyses still have large errors. To resolve this issue, a 
lot of theoretical ideas have been suggested. These ideas 
would be categorized broadly into two cases: modifica- 
tion of prior conventional models or claim of exotic states 
different from prior conventional structures. As one of 
improvements within the context of conventional quark 
models, it was shown that the DK mixing with the p- 
wave cs state plays an important role to lower the mass 
|l6j| . Apart from such a conventional picture, a possibility 
of the exotic configuration has been su gges ted: a four- 
quark state 0, |2£l , DK molecule |l7| , £ > s 7r atom 
[18f. and four-quark states with a DK mixing |21|. Us- 
ing perturbative QCD, Chen and Li calculated branching 
ratios of D*+ (2317) productions in the i3-meson decays 
and showed that the BELLE measurement of S-meson 
decays favors exotic multi-quark states of the D*+(2317), 
rather than conventional pictures |30j |. However, it was 
argued that four-quark states could not be bound for the 
charmed-strange mesons |3 1| . Thus, the theoretical ver- 
dict on this crucial issue is still not obtained. 

Another remarkable property of light-heavy systems is 
that in the heavy quark limit the j = 1/2 doublets at 
L = 1 could be interpreted as a chiral partner of those 
at L = [l|. According to this idea, one could 
expect that the mass-splittings between the + (l + ) and 
— (1~) states are almost equivalent to mj\r/3, where uin 
is a nucleon mass. The masses of the + and 1 + states 
obtained from this value are indeed near the values mea- 
sured by experiments. This property requires us compre- 
hensive discussion involving not only the excited states 
of + and 1 + but also the ground states of 0~ and 1~, 
when we investigate the masses of + and 1 + states in 
detail. 

Motivated by these, we formulate the QSRs and show 
the Borel-transformed forms for eight states of the cs and 
cn mesons with J PC = CP, 1~, + , 1 + in the next section. 
Note that we study these eight states independently in 
all the same framework and criteria for the BSR, and 
discuss their relative energy levels. 



III. FORMULATION OF BOREL SUM RULES 

We start by considering the covariant two-point 
function Tlpjg 2 ) with quark currents Jr(x) (R — 

o-i", o+i+) 

U R (q 2 ) = i [ d 4 x e ^(0|TJ fl (z)4(0)|0), (1) 



where |0) is nonperturbative QCD vacuum and T the 
time-ordered product. Jr(x) denotes the pseudo-scalar, 
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vector, scalar and axial-vector currents described respec- 
tively as 



the form of a step function: 



Jo- {%) 
J?_(x) 
J +(x) 
J? + (x) 



iq(x)j 5 c(x), 
q(x)y>*c(x), 
q(x)c(x), 
■q^ v q{x)^ v ^ b c{x), 



(2) 
(3) 
(4) 
(5) 



-Imn PH (s) 

7T 



F5(s 



-lmU OPE (s)9(s- s R ) 



(8) 



with q(x) being light-quark (d or u under the isospin sym- 
metry) fields or strange-quark (s) field and c(x) being 
charm-quark field at the point x, and rf ,v = q^q u /q 2 — 
g^ v . Note that, e.g., the Lorentz structure of the vector- 
current correlator with (x) and J^l (0) can be gener- 
ally expressed as the sum of the transverse invariant func- 
tion proportional to q^q v /q 2 — g^ v and the longitudinal 
invariant function proportional to q^q" /q 2 . We here con- 
fine our attention to only the transverse component. In 
the axial-vector case, automatically only the transverse 
component remains owing to addition of Lorentz tensor 
rf iV and hence we can get rid of the contribution from 
the pseudo-scalar meson in the lowest state |33| . 

We define invariant functions IT(q 2 ) (i = P,V,S,A), 
which are directly evaluated in our QSR, for respec- 
tive channels as follows: U -(q 2 ) = U P {q 2 ), U^(q 2 ) = 
(q^/q^-g^)U v (q 2 ), U 0+ (q 2 ) = n s (q 2 ) and 11^ (g 2 ) = 
[(fq v jq 2 — g^yU-Aiq 2 )- Under these definitions, the dis- 
persion relation satisfied by such invariant functions is 
given without any subtraction as 



U{q2) = 11 



ds 



Imn(s) 



(6) 



This equation relates the real part of the correlation 
function, which is valid in the high momentum region 
(— q 2 — > oo), with the imaginary part, which is valid in 
the low momentum region (q 2 — > Aq CD ) with a typical 
soft-scale Aqcd- Via this dispersion relation, we con- 
struct the QSR by equating the side of operator product 
expansion (OPE) and the phcnomcnological (PH) side 
related to the spectral function. The former is described 
as the product of Wilson coefficients and nonperturba- 
tive QCD vacuum condensates or quark masses, while 
the latter is parametrized by hadronic quantities such as 
resonance masses, couplings and the continuum thresh- 
old, etc. Thus, the QSR is represented in a simple form, 



ds W(s)- (imU 



PH, 



(s)-lmU OPE ( S )) = 0, (7) 



with the QCD continuum threshold Sr. F is a pole 
residue, where we take F — fRiri R k with Jr the cou- 
plings of the lowest resonances with respective parities to 
the hadronic current and tur a pole mass. The power k 
of m 2 R in the pole residue is taken to match the maximum 
power of s in the asymptotic s-behavior of the spectral 
function. For s > sr, we assume that the hadronic con- 
tinuum reduces to the same form with that obtained by 
an analytic continuation of the OPE, i.e. the perturba- 
tive terms, based on a hypothesis of the quark-hadron 
duality. Such an assumption is expected to smear the 
contributions of the higher radial excitations. This sim- 
ple ansatz would indeed reproduce well the characteristic 
features of the real spectrum in most applications of sum 
rules to mesons and baryons [32j| . 

Next, we display analytic forms of the QSR in each 
channel by implementing the Borel transformation to the 
QSR as in Eq. J7J 34J . Performing the OPE at dimension 
d < 6 operators, we obtain the relations of the BSR for 
four channels (0 + , _ , 1 + , 1~) of the cs meson: in the 
scalar + (pseudoscalar CP) channels, those read 
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(9) 



where W{s) is an arbitrary weight function, but must 
be analytic except for the positive real axis. The lower 
limit of the integration is given by the charm-quark Mass 

squared, m 2 . We adopt W(s) — exp(— s/M 2 ) in the BSR where g denotes the strong coupling, G 2 — G^G^ with 



below, where M is known as an artificial so-called Borel 
mass. We take the most general ansatz for the spectral 
function of the PH side, i.e., it would be saturated by one 
resonance in the narrow width limit and a continuum in 



gluon field G M „ and a ■ G = a^G^ . (O)o with the 
composite operator O represents (0|O(0)|0) with the local 
composite operator O(0) at the origin. On the other 
hand, in the axial-vector 1 + (vector 1~) channels, we 
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similarly obtain 

^2 2 -m 2 ,/M 2 

/ 1± m 1± e i±' 

8tt 2 7 m 2 V s / V s 

3m c m s s 4a s (s) 



it 



-m 2 a /M 



(2s + ml){s - ml) 3 tt 
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1 TfL^ 

T --^|m c (s 5 er- Gs) 



1/^G 2 

12 \ 7T / 

11/ m 2 mf. , 



207r 1 / m 2 1 m" L 



(10) 



Here, we approximate at the first order of m s , because 
m s is small enough compared to m c or a typical scale 
of M (M ~ 1 GeV). The 0(a s ) corrections to the per- 
turbative contributions are given as the function Rq(x) 
for the (pseudo-) scalar channel |35l l36f and the function 
Ri{x) for the (axial-) vector channel [36ll37| as follows: 



9 3 1 — x 

Ro(x) = - + 2Li 2 (x) + lnxln(l - x) - -In 

4 2 x 



— ln(l + i In 



l-x 



1 



-lnx 



(11) 



13 3 1 — x 

R x {x) = h 2Li 2 (x) + lux ln(l — x) In 

4 2 x 



ln(l — a;) + x In 



1-z 



1 -x 



-lnx 



(3 + i)(1-i), l-x 
m 



2x 



2 + x 
5 2x 



x (2 + x){l-x) 2 
2x 



2 + x 2 + x (2 + x)(l-x) 



lnx 
(12) 



with the Spence function Li 2 (x) = — Jjf cfti _1 ln(l — t). 
As for the running coupling constant a s {s) appearing in 
the perturbative terms of Eqs. 10 and l]lUp. we approx- 
imate it by a one-loop form, a s (s) — 47r/(91n(s/AQ C£) )) 
with Aq C£) = (0.25 GeV) 2 , which is determined to repro- 
duce a s (lGeV) ~ 0.5 

As easily seen from Eqs. © and (|10fl . we can get the 
relations of the OPE in the 0~ and 1 _ channels by re- 
placing simply m c by —m c for those in the + and 1 + . 
Also, replacing m s and ss by m n and nn respectively, we 
can derive the corresponding BSRs for the non-strange 
cn mesons except for the only difference of the continuum 
thresholds. We here take the massless limit (m n = 0) of 
light quarks for the BSR of the cn mesons. These re- 
sults indeed reproduce the corresponding BSRs for light 
mesons in the vector, scalar and axial-vector channels 
(sec the results of Rcf . 39] in vacuum) . Looking on the 
relations of OPE in the right-hand sides of © and (II Ul) . 
we find that the difference of property between the cn 



and cs mesons would be caused by the effect of strange 
quarks that is mainly the difference between lowest di- 
mensional condensates, m c (nn) and m c (ss) . For the 
pseudoscalar channel of the cn mesons, the similar rela- 
tion was used in Ref. |40j, although their relations are 
somewhat different from ours for whole sign of a term 
including gluon condensate and a part of terms includ- 
ing four-quark condensates. This difference, however, is 
insensitive to the final results. For the axial and vector 
channels of the cn mesons, under approximation neglect- 
ing gluon condensates and four-quark condensates, the 
same relation was used in Ref. [41| . More recently, Nari- 
son presented explicit BSR forms in the scalar and pseu- 
doscalar channels of the cn and cs mesons 23] . They are 
very similar to our forms. 

In order to extract the values of physical quantities 
from the BSR, we use the following standard values for 
QCD parameters appearing in the OPE: m s — 0.12 GeV 
[Hi, m c = 1.46 GeV & (nn)o = (-0.225 GeV) 3 [H, 
(ss) = 0.8x(nn) GeV 3 !^, <^G 2 ) = (0.33GeV) 4 x2 
HI, (figa ■ Gn) = M 2 (nn) with M 2 = 0.8 GeV 2 [H, 
Gs) n = M 2 (ss) , a s (nnf = 0.162 x 10~ 3 GeV 6 



[sga 



and i 



= 0.8 2 



GeV 6 . Note that we eval- 



) = U.S" x a s {nn) n 
uated four-quark condensates in terms of (nn) Q , suppos- 
ing that the factorization or vacuum saturation in an in- 
termediate state of four-quark matrix element can work 
well in a good approximation |22j |. Although the ambi- 
guity of the absolute value of the four-quark condensate 
still remains large as discussed by several authors (see 
Refs. [3(1 E3), the detailed discussion is not much rel- 
evant to our present analysis. It is because the contri- 
bution from higher dimensional condensates is negligible 
compared to the lowest dimensional operators (d = 3) 
and thus the final results are independent of details of 
four-quark condensates. 

Finally, taking a logarithmic derivative with respect 
to 1/M 2 for both sides of © and (JTDJl, we get rid of 
/^-dependence completely from their equations and thus 
can derive the Borel curve of the resonance masses m 2 R 
as functions of M 2 and sr. 



IV. NUMERICAL RESULTS 

We determine plausible values of the resonance mass 
extracted from the Borel curves in accordance with the 
following criteria. Its method is in order: 

1. We plot the Borel curves of the resonance mass niR 
as a function of M with taking some selected val- 
ues of the continuum thresholds ^/sr. In general, 
the values of the thresholds might be different from 
the phenomenological value of the first radial exci- 
tation mass and rather would be smaller. This is 
natural, because our simple ansatz in the PH side of 
the QSR cannot take into account the complicated 
structure of the resonances in this region. We, 
therefore, set the maximum value of ^/sr, ^/s™ x , 
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to be the mass of the first radial excitation pre- 
dicted by the potential model 

2. In the OPE side, the contribution of the nonpertur- 
bative corrections should be generally small com- 
pared to that of the perturbative terms in large 
momentum region. Since we, however, truncate 
the OPE at d = 6 by hand, we would make the 
impact of this approximation as small as possible. 
In naive sense it is satisfied requiring that the con- 
tribution from higher dimensional operators than 
d — 6 becomes thoroughly smaller than the pertur- 
bative terms in large region of M . In such a region 
of M, the convergence of the OPE can indeed be 
ensured to some extent. Here, we quantitatively 
require a constraint that the effect of d = 6 term 
is less than 20% compared to that of the leading 
terms in Eqs. J^J and QlOfl. This assures us to get 
physical quantities with about 10 ~ 20% accuracy. 
This condition gives the lower limit of M, M m i n . 

3. In the PH side, we would make the continuum con- 
tribution as small as possible, since we are inter- 
ested only in the pole mass. It would be possible 
by taking small region of M. We require that the 
continuum contribution is less than 20% compared 
to that of the perturbative contribution. For the 
leading terms in Eqs. © and IjlOfl . the continuum 
contribution corresponds to the integration over s 
from sr to infinity, while the perturbative contri- 
bution has the integration from m 2 c to infinity. This 
constraint also assures us to get physical quantities 
with about 10 ~ 20% accuracy. From this con- 
dition, we can give the upper limit of M, M max , 
which strongly depends on the threshold in con- 
trast to M m i n . 

4. Thus, we set the Borel window, M m i n < M < 
M max , so that the BSR would work well within 
the window and we extract more reliable values of 
physical quantities from the Borel curves. 

5. For sfs~R < \Js^ ax , we look for the stable regions 
(Borel stability) of the Borel curves vs. M within 
the Borel windows. Then, from the curve having 
the widest Borel stability, we read off the resonance 
mass and the threshold as the most reliable values 
of physical quantities. In general, this would give 
the central values of physical quantities under con- 
sideration. Finally, we adopt all resonance masses 
and thresholds to the extent that the Borel stabil- 
ities on the curves exist within the Borel windows. 
This would give errors for the central values. Note 
that in actual evaluation for the Borel stability we 
do not mind its accurate positions and widths of 
the stable regions. Such an ambiguity also should 
be included in the above errors. 

Following the above criteria, we plot the Borel curves 
of resonance masses iur vs. M in two ground states 




0.8 0.9 1.0 1.1 1.2 1.3 1.4 

M [GeV] 



FIG. 2: The Borel curves of m D ^ Q -^ vs. M, where ^/s.d(o-) = 
2.3 ~ 2.6 GeV at 0.05 GeV intervals. Explanation for the lines 
is given in the text. 

(R = 0~,1~) of the cn mesons. We first study these 
two channels to determine a size of charm-quark mass 
m c , which is treated as an input data in the OPE side. 
Since the quark mass is a scale-dependent quantity, its 
size runs with the energy scale of the system, e.g. we 
adopt m c = 1.2 ~ 1.5 GeV in the QSR approach [44Ll45| . 
We determine m c to reproduce the masses of both states 
as well as possible. This leads to m c = 1.46 GeV as listed 
before. Using this input, we shows the curves for the _ 



2. 1 2 1 1 1 1 1 1 i i i 




1.0 1.1 1.2 1.3 1.4 
M [GeV] 



FIG. 3: The Borel curves of m^^-j vs. M, where ^s D ^- ) = 
2.45 ~ 2.65 GeV at 0.05 GeV intervals. 

channel with y ^Fj = 2.3 - 2.6 GeV at 0.05 GeV in- 
tervals in Fig. [21 and for the 1~ channel with ^/socx-) = 
2.45 - 2.65 GeV at 0.05 GeV intervals in Fig. E3 In 
Fig. |2 left long (right short) vertical-lines on the curves 
represents the lower (upper) limit of the Borel window, 
M m i n (M max ) , which very weakly decreases (strongly in- 
creases) with the threshold. We show the Borel stabilities 
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2.55 



2.50 



2.45 



2.40 



FIG. 4: The Borel curves of m 




TABLE I: Numerical results of the resonance mass run and 
the continuum threshold y/~SR. They are listed for four chan- 
nels (R = 0~,1~,0 + ,1 + ) of the cn mesons. We refer to 
Ref. |26 
from Ref. 9 



for the experimental data and the predictions 
for maximum values of the thresholds. 



R 


m R 


D [GeV] 
m fl (exp.) 




pt. model \9j 


Q- 


1.90 ± 0.03 


1.869 


2.45 ± 0.15 


2.589 


1" 


2.00 ± 0.02 


2.010 


2.55 ± 0.05 


2.692 


0+ 


2.45 ± 0.03 


2.351 


2.90 ± 0.10 


2.949 


1+ 


2.38 ± 0.05 


2.438 


2.70 ± 0.10 


3.045 



m D{1+) = 2.38 ± 0.05 GeV, which are read off Figs. H 
13 respectively. Here, we took ^/s£>(o+) = 2.8 ~ 3.0 GeV 



D(0+) 



M 



2.8 ~ 3.0 GeV at 0.05 GeV intervals. The lower limits of the 
Borel windows are below M = 1.1 GeV. 



by dotted lines within these Borel windows. The length 
of the lines corresponds to the range of stable regions. In 



where y ^Dfo+j = anc ^ y/ s D(i+) — 2-5 ~ 2.8 GeV, respectively. We summa- 




FIG. 5: The Borel curves of m D ^+^ vs. M, where ^Sd(i+) = 
2.5 ~ 2.8 GeV at 0.1 GeV intervals. 



the 0~ channel, we cannot find any Borel stability below 
y /s D(o-) = 2.3 GeV, while above the value there exist 
such stable regions on the curves until the upper limit 
y/ s D{o-) — 2.6 GeV. Thus, we get the resonance mass, 
m D(o-) — 1-90 ± 0.03 GeV. In the l - channel, we can- 
not find any Borel stability below ^/sd(i-) = 2.45 GeV, 
while above the value there exist such stable regions on 
the curves until ^/sn(i-) — 2.6 GeV, which is below the 
upper limit ^Jsd\i-) — 2.7 GeV. Thus, we get the res- 
onance mass, m D (i-} = 2.00 ± 0.02 GeV. Indeed, the 
obtained results for resonance masses are near experi- 

1.867 GeV and m^Z. ~ 2.010 



mental values, "i^,'-) 

GeV, respectively. In a similar analysis, we get the re- 
sults of other channels, mr>(o+) = 2.45 ± 0.03 GeV and 



rize these results in Table H] From the table, we find that 
the excited (p-wave) 1 + state reproduce the experimental 
data within errors, while the excited + state seems to be 
overestimated by about 100 MeV in comparison with the 
experimental data. However, if we take into account the 
broad width (Tjj^Q+j ~ 329 MeV) measured by experi- 
ments, our result could be not in conflict with such data 
(within the large width). Also, the obtained thresholds 
are different among the four channels and thus have a 
level structure similar to the obtained resonance masses; 
the order of this energy level is the same as that of the 
resonance masses. 



2.05 

q 

# 2.00 




0.8 0.9 1.0 1.1 1.2 1.3 
M [GeV] 



FIG. 6: The Borel curves of m Dt 



(o-) 



vs. 



M, where 



2.3 ~ 2.7 GeV at 0.05 GeV intervals. 



Next, we discuss the case of the cs mesons following 
the above criteria. Here we fix m s = 0.12 GeV as listed 
before. Indeed, final results of the resonance masses are 
almost insensitive to the adopted values of m s . For exam- 
ple, the result with m s = 0.15 GeV changes by less than 
0.5 % of that with m s = 0.12 GeV. Hence, the error com- 
ing from such a mass ambiguity is smaller than the error 
originated from the criteria mentioned above. We show 
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0.9 1.0 1.1 1.2 1.3 1.4 1.5 
M [GeV] 



FIG. 7: The Borel curves of m 



V S -D S (1-) 



D,(l-) vs - 

= 2.4 ~ 2.9 GeV at 0.1 GeV intervals. 



M, where 



give the mass-splitting of the resonance masses. We 



2.50 




the results of the channel in Fig. The behaviors of 
the curves and the Borel windows are very close to those 
of the D{0~) in Fig. We get m^o-) = 1-94 ± 0.03 
GeV for ^/Sd 3 (o-) = 2.3 ~ 2.7 GeV by applying the 
above criteria. This value is almost consistent with the 
experimental data (1.969 GeV) within errors and higher 
by 40 MeV than our result for the cn meson. Thus, our 
mass-splitting between the D s (0~) and -D(0~) is underes- 
timated by about 100 MeV in comparison with the mea- 
sured one, which would just correspond to a deviation 
between light-quark and strange-quark masses. Another 




x 2.48 



1.3 1.4 1.5 
M [GeV] 

FIG. 8: The Borel curves of m Dg ^ 0+ ^ vs. M, where 
^ s ga(0+) = 2.85 ~ 3.05 GeV at 0.05 GeV intervals. The 
lower limits of the Borel windows are below M = 1.15 GeV. 

remarkable point is that there exists also a deviation be- 
tween their thresholds as well as the mass-splitting of 
the resonance masses. This deviation is near the value 
of the mass-splitting. It is, therefore, expected that such 
a deviation for the thresholds plays an essential role to 



FIG. 9: The Borel 

V S D S (1+) = 2 - 7 



curves of m,o s ci_+) vs. M, where 
2.95 GeV at 0.05 GeV intervals. 



also show numerical results of other three channels in 
Figs. 0|HllS|in the same analysis. Every graph has the 
curves similar to the corresponding channel of the cn me- 
son, respectively. In the 1~ channel, the obtained mass 
m Ds (i-) = 2.05 ± 0.04 GeV for y ^s^F) = 2.4 - 2.9 
GeV is somewhat smaller than the experimental data 
(2.112 GeV) like the case of the cn meson. In this channel 
also, our mass-splitting between the D a (l~) and D(l~) 
is about a half of the experimental observation, i.e. 
~ 50 MeV, and also a difference of the obtained thresh- 
olds between their mesons is non-zero, i.e. about 100 
MeV. In the + channel, we find that the obtained mass 
m Ds{Q+) = 2.48 ± 0.03 GeV for ^ s Ds{0+) = 2.85 - 3.05 
GeV is overestimated by about 160 MeV compared with 
the experimental data (2.317 GeV), as seen in the case of 
the D(0 + ). This result is very close to that from the po- 
tential models H,E1- Contrary to the case of the D(0 + ), 
the D s (0 + ) meson has been observed as the very nar- 
row state (r £ ) s ( +) < 4.6 MeV). Therefore, this deviation 
(160 MeV) is significant and our result supports conclu- 
sions from Ref. H, 0, 0, 0, . A deviation between 
the thresholds of the D s (0+) and D(0+) is about 100 
MeV. In contradiction to the result in Ref. 0, |t| , how- 
ever, our result in the 1 + channel, m-D s (i+) = 2.41 ±0.05 
GeV for ^sn,(i+) = 2.7 ~ 2.95 GeV, is in agreement 
with the experimental data (2.459 GeV) within errors. 
Such agreement with the data could be due to disregard 
of mixing effects between J p (°' = ]+(+) and al- 
though Refs. 0, have taken into account the effect. 
Here, our result for the threshold gives almost the same 
value with the case of D(l + ). These final results for four 
channels of the as mesons are summarized in Table ITT1 
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TABLE II: Numerical results of the resonance mass run and 
the continuum threshold %/sr. They are listed for four chan- 
nels (R = 0~,1~,0 + ,1 + ) of the cs mesons. For compari- 
son, we attach experimental average values observed recently 
Elk&QllMbl m each channel and the masses of the first radial 
excitations, which is predicted in Ref. 0. 







D s [GeV] 






R 


ma 


rafl(exp.) 


\fsa 


pt. model \9] 


(T 


1.94 ± 0.03 


1.969 


2.50 ± 0.20 


2.700 


1" 


2.05 ± 0.04 


2.112 


2.65 ± 0.15 


2.806 


0+ 


2.48 ± 0.03 


2.317 


3.00 ± 0.15 


3.067 


1+ 


2.41 ± 0.05 


2.459 


2.70 ± 0.25 


3.165 



V. SUMMARY AND DISCUSSION 

We have presented detailed calculations of the masses 
for the 0~,1 _ ,0 + ,1 + charmed-mesons (D, D s ), moti- 
vated by recent discovery of very narrow two states 
(£1^(2317), L»^(2459)) by the BABAR, CLEO, BELLE 
and FOCUS Collaborations. Applying the original QSR 
techniques of SVZ [22J for the conventional light-heavy 
systems, we have performed the OPE up to d = 6, in- 
volving corrections to the order a s and to the order m s 
(or m„). In this method we have taken into account full 
order of the l/m c -expansion, because m c is not a large 
mass enough to approximate by the first few terms of 
the expansion. We have analyzed their masses using the 
Borel-transformed QSR, because this analysis can care- 
fully deal with the continuum contribution or the thresh- 
old dependence, if the resonance mass, which we would 
derive from this analysis, is sensitive to such a contri- 
bution Especially, such care is important for dis- 
cussion of the spectroscopy associated with (hyper-)fine 
structures as in our considering four channels. In fact, as 
seen in Tables HI and HTl we found that the mass-splittings 
among four channels of the cn or cs mesons and those in 
the same channels between the cn and cs mesons can be 
largely affected by the values of their thresholds. In or- 
der to derive a reliable resonance mass from the threshold 
dependence, we performed numerical calculations by im- 
posing stringent criteria listed in sec. IV for the BSR. We 
compare our final results for resonance masses of cfi and 
cs mesons with those of the experimental data in Fig. ITU1 
From this figure, we find that the mass of the D s (0 + ) 
meson is overestimated by about 160 MeV in compari- 
son with the experimental data, which is consistent with 
the potential model analysis [3, @ • Although a similar 
tendency is also seen in the mass of the cn meson in 
the same channel, these results could be consistent with 
the experimental data within the measured broad width. 
On the other hand, such a tendency is not seen in other 
channels. Conversely in other channels the cs masses are 
underestimated somewhat in comparison with data, in- 
dependent of strange-quark mass adopted in our calcula- 
tions. Hence, the 1 + — 1~ and the + — mass-splittings 
are unequal in disagreement with the chiral multiplet 
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- 2459 

2438 

2351 
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™ ~01869 
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• O Data 
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n-(+) ,"(-) „+(+) ,+(+) T P(C) 
"l/2 A l/2 "l/2 A l/2 J j 

1S 3P„ SPj 2 S+ l Lj 

FIG. 10: Spectroscopy of the cn and cs mesons with J p — 
0~,1~,0 + ,1 + . The solid lines for the cs mesons and dashed 
lines for the cn mesons show our numerical results. The cor- 
responding experimental data are given by closed and open 
circles, respectively, as in Fig. Q 



structure predicted by chiral effective theory for heavy 
mesons. It might be expected that the measured low 
mass of D* + (2317) is a manifestation of any exotic states 
with the structure of a four-quarks or a molecule. In 
fact, we have seen very recently |2|| that the assignment 
of the -0*^(2317) to the charmed-strange four-quark me- 
son, F\, with {1,13) = (1,0) is favored by the ratio of 
the measured rates for the D*+(2317) -> to the 

D s *+(2317) - Dfir° decay i, T(D*+7)/T(D+tt ) < 
0.052. It is, however, hard to reconcile its assignment 
to an isosinglet state (the conventional {cs} or an isos- 
inglet four-quark or a molecule) with the experimental 
constraint. Our conclusion for the charmed mesons with 
+ is different from that of Ref. |23, where the similar 
QSR analysis was performed only for 0~ and + states 
of the cs and cn-mesons without using stringent criteria 
like in our analysis. 

In this work, we neglected renormalization-group im- 
provement for the currents and the operators in the con- 
densates, which could be important to improve the Q 2 
range of the validity of the derived QSR. Also, as is well 
known, the physical quantities extracted from the QSR 
analysis generally have theoretical ambiguity of about 
10 ~ 20 %, which comes from the condensates and our 
criteria, etc. However, these quantitative ambiguities 
would not qualitatively change our conclusion as men- 
tioned above, because we have evaluated totally eight 
states of the cn and cn mesons within the same method 
and criteria. In such an analysis we found an anomalous 
feature for only the charmed-strange scalar meson. 
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